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ABSTRACT 



Nuclear obscuration plays a key role in the initial phases of AGN growth, yet not many highly obscured AGN are currently known 
beyond the local Universe, and their search is an active topic of research. The XMM-Newton survey in the Chandra Deep Field 
South (XMM-CDFS) aims at detecting and studying the spectral properties of a significant number of obscured and Compton-thick 
(A^H ^ 10^"^ cm~^) AGN. The large eff'ective area of XMM-Newton in the 2-10 and 5-10 keV bands, coupled with a 3.45 Ms nominal 
exposure time (2.82 and 2.45 Ms after lightcurve cleaning for MOS and PN respectively), allows us to build clean samples in both 
bands, and makes the XMM-CDFS the deepest XMM-Newton survey currently published in the 5-10 keV band. The large multi- 
wavelength and spectroscopic coverage of the CDFS area allows for an immediate and abundant scientific return. In this paper, we 
present the data reduction of the XMM-CDFS observations, the method for source detection in the 2-10 and 5-10 keV bands, and the 
resulting catalogues. A number of 339 and 137 sources are listed in the above bands with flux limits of 6.6 x 10~^^ and 9.5 x 10~^^ erg 
s~^ cm~^, respectively. The flux limits at 50% of the maximum sky coverage are 1.8 x 10~^^ and 4.0 x 10~^^ erg s~^ cm~^, respectively. 
The catalogues have been cross-correlated with the Chandra ones: 315 and 130 identifications have been found with a likelihood- 
ratio method, respectively. A number of 15 new sources, previously undetected by Chandra, is found; 5 of them lie in the 4 Ms area. 
Redshifts, either spectroscopic or photometric, are available for ~ 92% of the sources. The number counts in both bands are presented 
and compared to other works. The survey coverage has been calculated with the help of two extensive sets of simulations, one set per 
band. The simulations have been produced with a newly-developed simulator, written with the aim of the most careful reproduction 
of the background spatial properties. For this reason, we present a detailed decomposition of the XMM-Newton background into its 
components: cosmic, particle, and residual soft protons. The three components have diff'erent spatial distributions. The importance of 
these three components depends on the band and on the camera; the particle background is the most important one (80-90% of the 
background counts), followed by the soft protons (4-20%). 

Key words. Catalogs - surveys - galaxies: active - methods: data analysis - X-rays: general 
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1. Introduction 



^ The understanding that the X-ray background (XRB) is due to 
unresolved emission by discrete extragalactic sources ([Giacconi 



|& Zamorani|1987|[Schwartz et aLlf ^Tgi'Macca caro et al.|199T] )7 
and that obscured AGN are a key ing redient to produce the 
observed shape of th e XRB spectrum ( |Setti & Woltjer||1989| 
[Comastri et al.||1995| ) has put a strong focus on the role of ab- 
sorption as a driver of the observable properties . This idea has 
shaped th e original model of AGN unification (Anto nucci|1993[ 
|Urry & Padovani__1995} and keeps its importance unaltered in 

* Based on observations obtained with XMM-Newton, an ESA sci- 
ence mission with instruments and contributions directly funded by 
ESA Member States and NASA. 



the current modelling of AGN and of their evolving populations 



( Treister & Urry|2Q05[prii et al.|2007[|Treister et al.|20Q9a| ) 



Nuclear obscuration might be associated with the initial 
phases o f AGN em ission ( [Page et al.|2QQ4t [Hopkins et al.|2Q"Q6l 
[Menci et al.|20Q8| ). The idea is that a large gas reservoir is avail- 
able at high redshift to feed (and obscure) an accreting super- 
massive black hole. The same gas reservoir would also sustain 
high star formation rates in the host galaxy, leading to the broad 
similarity between the cosmic histories of accretion and star for- 
mation. 

Deep X-ray surveys are a primary tool for the census of AGN 
and the study of the properties of their populations. Since the 
launch of the Einstein observatory, the first imaging X-ray tele- 
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scope, many surveys have been perfomed with different combi- 
nations of area and flux limit and in many locations on the sky. 

The Chandra Deep Field South (CDFS) survey started with 
1 Ms of observatio ns ([Giaccon i et al. 2002 ) and was afterwards 



20 



2001-2002 



2008-2010 



extended to 2 Ms ( |Luo et aLj 2008 ) and finally to 4 Ms ( |Xue 
et al.| 2011 hereafter Xll). It currently reaches the deepest X- 



-1 



ray fluxes ever probed, with a flux limit of 9.1 x 10" erg s 
cm"^ in the soft (0.5-2 keV) band. However, given the energy 
dependence of its efl'ective area, Chandra is much less sensitive 
in the hard and very hard bands (2-10 keV and 5-10 keV, re- 
spectively) which are the most important for obscured objects. 
The flux limit for 2-10 keV is currently 5.5 x 10"^^ erg s"^ cm"^ 
(Xll). The only analysis for the 5-10 keV band was done for 
the initial 1 Ms survey ( Rosati et al.||2002 ) and yielded a flux 
limit of 1.2 X 10"^^ erg s^ cm"^; however, the 4 Ms survey was 
analysed in th e 4-8 keV band, mo ving the flux limit to 6 x 10"^^ 
erg s"^ cm"^ ( Lehmer et al. 2012| ). A larger area centred around 
the CDFS was surveyed with Chandra and with much shorter 



exposure times as the Extended CDFS (ECDFS; |Lehmer et al. 
j2005j). 

XMM-Newton is better suited to collect photons in the 
hard and very hard bands, because at energies > 5 keV the 
XMM-Newton efl'ective area drops less sharply than Chandra'^. 
However, it has a larger point spread function (PSF), which 
makes it sufl'er more from source confusion, and it has a higher 
instrumental background. While XMM-A^^w^6>^ surveys cannot 
probe very faint fluxes as Chandra, they provide good quality 
spectroscopy in the whole 0.5-8 keV band for a large num- 
ber of sources. The large efl'ective area of XMM-Newton has 
been exploited in several surveys, including the Lockman Hole 
(with 1 Ms of exposu re, Brunner et al. 2008), the ELAIS-Sl field 
pQg]), COSMOS ( ICappelluti et al.|2007||2QQ9l ), 



(Puccetti et al. 



and 2XMM ( jMateos et al.|2008] ). 



The XMM-Newton survey in the CDFS (hereafter XMM- 
CDFS) has been started with the main aim of finding obscured 
AGN, especially the Compton-thick ones, and studying their 
properties. With a nominal exposure of 3.45 Ms centred on a 
single point on the sky it is currently the deepest XMM-Newton 
observation ever performed. Moreover, it has the advantage of 
having almost pan-chromatic coverage and a large number of 
available optical spectra and redshifts. 

The first results of the XMM-CDFS survey have been pub- 
P011| ). Diff'erent approaches to the 



lished in Comastri et al. 



selection and study of obscured sources are being e mployed 
( Iwasawa et al.|[2QT2{ |Georgantopoulos et al. |2013[ spectral 
stacking (Falocco et al., submitted to A&A) and spectral anal- 
ysis of a flux limited sample (Comastri et al., in prep.) are also 
ongoing. Extended sources will be discussed in Finoguenov et 
al. (in prep.). 

In the present paper, we present the hard and very hard cata- 
logues, and the Log A^-Log S in the two bands. A number of 339 
and 137 sources are detected, respectively, with a significance 
threshold roughly equivalent to 4.8cr. Supplementary catalogues 
are provided for sources with lower significance (4(t). 

In Sect. |2] we present the observation details and the data re- 
duction procedure. In Sect. [3] we describe the source detection 
procedure and present the catalogues in the hard and very hard 
bands. The survey coverage and number counts are presented in 
Sect.m Simulations of mock XMM-CDFS fields are instrumen- 
tal in determining the coverage, and are presented in Sect. [5] 
The simulations also used in Sect. |6] for a first estimate of the 
number of spurious sources. We present in Sect. |7] an analysis 
of the source confusion in the hard band. In Sect. [8] we iden- 
tify the Chandra counterparts to the XMM-Newton sources. We 
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Fig. 1. Total surface brightness in the 2-10 keV band for the 33 
obsids of the XMM-CDFS, after excluding th epositions around 
the 30 brightest sources. The increase in background intensity 
between 2002 and 2008 is evident. The obsids are numbered 
from 1 to 33 according to observation date; the vertical dotted 
line marks the separation between the 2001-2002 and 2008- 
2010 observations. Filled triangles: PN; filled circles: MOSl. 
The M0S2 camera has values very similar to the MOSl. 

summarize our conclusions in Sect. [9] Finally, in Appendix |A| 
we present the simulator that we have developed. 

2. Observations and data reduction 

The Chandra Deep Field South was initially observed by XMM- 
Newton in the years 2001-2002 (P.I.: J. Bergeron) with an expo- 
sure of 541 ks. It was observed again, for the proposal which 
led to this series of papers, in the years 2008-2010 (P.I.: A. 
Comastri), to reach a nominal exposure of 3.45 Ms. 

Table [T] shows the observation dentification number (here- 
after obsid), date, pointing and exposure for all observations. 
The number of obsids is 33; considering the MOSl, M0S2 and 
PN cameras, the XMM-CDFS survey comprises a total of 99 
event files. We used the XMM-Newton SAS software for our 
analysis (version 10 for the initial processing and the catalogue; 
version 11 for the simulations; there have not been significant 
changes between the two versions which could aff'ect this work), 
with the help of many custom-developed scripts to automate 
most of the data processing. 

2.1. Background flares and quiescent level 

Each one of the 99 event files was screened individually for 
background flares. Full-field light curves were generated in the 
10-13.5 keV interval. A 3cr-clipping procedure was applied to 
the light curves to identify and reject the high-background pe- 
riods. This procedure worked well most times, though it failed 
for some event files which were severely flared. In such cases, 
we adopted the nominal count rate thresholds given by the 
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Fig. 2. Images of the XMM-CDFS in the 2-10 (left) and 5-10 keV (right) bands. The colour wedges show the total counts as the 
sum of data from the MOSl, M0S2 and PN cameras. 



XMM-Newton documentatior[^ to identify and reject the high- 
background periods; these observations are marked in Table [l] 

A relevant and unexpected feature is an increase by a fac- 
tor of ~ 2 in the quiescent background level in the 2008-2010 
observations with respect to the years 2001-2002 (Fig.[T]). The 
increase is seen in all three cameras and is due to the instrumen- 
tal ("particle") component of the background (see Sect. [5]). The 
reason for this increase is not clear, though one possibility is that 
it is related to the Solar cycle. 

The images in the 2-10 and 5-10 keV bands, showing the 
total counts from the MOS and PN cameras, are presented in 
Fig. [2] The total exposure, after the high-background period fil- 
tering, is thus 2.82 and 2.45 Ms for the two MOS and the PN, re- 
spectively. The exposure maps for the 2-10 and 5-10 keV bands 
are shown in Fig. [3] 

2.2. Copper lines complex 

For the PN camera, the strongest background feature in the 2-10 
keV interval is a complex of lines (Cu Ka, Ni Ka, Zn Ka, the 
Cu being the strongest one) around 8 keV, due to fluorescence 
from the mirror holding structure, and which alone can make 
~ 30% of the total counts in the 2-10 keV band. This complex 
has a well-defined spatial pattern with a central region where it 
is virtually absenj^ 

Throughout this paper, in PN data we have excluded the en- 
ergy intervals pertaining to this complex. We have chosen dif- 
ferent trade-off's for the 2-10 and 5-10 keV bands. For the 2-10 
keV band, we have excluded the Cu complex (energy intervals 
7.2-7.6 and 7.8-8.2 keV) across all the field, to favour a uniform 



^ http://xmm.esac.esa.int/sas/current/documentation/threads/ 
EPIC_filterbackground. shtml 

2 XMM-Newton Users Handbook, Issue 2.10, Sect. 3.3.7.2; 
[http://xmm.esac.esa.int/ external/ xmm_user_support/ documentation/ 
uhb/ epicintbkgd.html 



coverage. For the 5-10 keV band, given the smaller energy band 
and the lower signal/noise ratio, we have been more aggressive 
in excluding the Cu complex (energy intervals 7.35-7.60, 7.84- 
8.28 and 8.54-9.00 keV), but we have done so only in the outer 
detector regions where the line is present. Our aim for 5-10 keV 
has been to retain the photons around 8 keV in the centre of the 
FOV at the expense of a less uniform coverage of the field. 



2.3. Pointings and astrometry 

The XMM-CDFS observations have been taken with some dif- 
ferences in the pointing direction, in order to have a field cov- 
erage as uniform as possible with the PN detector. Because 
of the field visibility, the observations have been performed in 
the months of January, February and July. The orbit of XMM- 
Newton modulates the position angle, so that the satellite ro- 
tates along its optical axis by about 180 deg every six months. 
Therefore, when looked at in detector coordinates, observations 
performed in summer appear "upside-down" with respect to the 
winter ones. While the MOS detectors are centred on the respec- 
tive telescope optical axes, the PN is shifted by 2\ The pointing 
coordinates have thus to be shifted by twice that amount in or- 
der to position the PN on the same sky region in the two seasons. 
This, however, has the eff'ect of making the MOS detectors offset 
with respect to the nominal CDFS coordinates. 

Minor shifts of < V were added to all observations, to 
smooth the chip gaps and to avoid the possibility that any source 
could fall in a gap for a major part of the survey. 

On top of these patterns, there are smaller random errors 
(~ V\ with 4 obsids in which they were ~ 4'') in the point- 
ing direction. These have been estimated by running a source 
detection on the inner 5' of each obsid, and cross-correlating the 
results with the Chandra catalogue (XI 1) using a modified ver- 
sion of Chandra's routine align_evt. The event and attitude 
files have been corrected for these shifts. All coordinates have 
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Fig. 3. Exposure maps of the XMM-CDFS in the 2-10 (left) and 5-10 keV (right) bands. The Chandra 4 Ms and ECDS areas, and 
#ie pointing positions and groups are superimposed with different symbols and letter marks according to the observing season (see 
Table[T]). The colour wedges show the exposure time in s as the sum of the MOSl, M0S2 and PN exposures. 
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been aligned to those from XI 1 (which was itself aligned to VLA 
radio frame). 

The XMM-Newton PSF has a ~ 5'' FWHM on-axis, which 
degrades at large off-axis angleq^ We measured the observed 
2-10 keV FWHM of the 50 brightest sources for the total 3Ms 
image and we found that the median FWHM, across the whole 
field, is 8.5''. 

2.4. Colour image 

A colour image of the XMM-CDFS has been derived from the 
cleaned, astrometry-aligned data. Three images in the 0.4-1, 1-2 
and 2-8 keV energy intervals have been produced. To better en- 
hance the sources, the background has been suppressed up to dif- 
ferent thresholds in the three bands. Finally, the non-linear scal- 
ing and colour-mixing procedure by Luptonetal. (2004), whose 
main advantage is that bright sources are not "whited-out", has 
been applied producing the colour image shown in Fig.|4] 



3. Catalogue 

The 2-10 and 5-10 keV catalogues have been produced inde- 
pendently. A two-stage process has been used: the first step 
is a detection with a wavelet code (PWXDetect), followed by 
maximum likelihood fitting (EMLDetect) of the sources. This 
detection method can be regarded as a variant of the standard 
XMM-Newton detection process, which uses a sliding box de- 
tection algorithm for the first stage instead of the wavelet (see 
e.g. Brunner et al.|[2008, and |Cappelluti et al.||2009 for applica- 
tions to the Lockman Hole and COSMOS surveys, respectively). 
The advantage of using a wavelet-based tool is that it produces 
more accurate source positions (Pucce ttTet al.|f2009) , and that 
it is especially suited for crowded fields, as happens in a few 
small areas of the CDFS. Another diff'erence with respect to the 
Lockman Hole and COSMOS applications is that we chose a 
higher threshold (i.e., more conservative) for the first stage than 
for the second (see below); doing so means that the selection is 
mostly done at the first stage, and that the second is used mainly 
for the counts estimate. 

First, we ran the Palermo Wavelet XMM Detect (PWXDetect) 
software to identify so urce candidates. P WXDetect is a version 
of the PWDetect tool ( [Damiani et al.|19 97), specially tuned for 
XMM-Newton data. It can combine data from the three EPIC de- 
tectors and from diff'erent observations, even taken with diff'erent 
pointings, provided that the pointing boresights are close enough 
that a single model for the off'-axis angle dependence of the PSF 
can be used (in practice, this amounts to a few arcmins). We used 
this feature to perform the detection on all obsids and cameras. 
We chose a significance threshold for source detection at the 4(T 
level, corresponding to a probability of 6.3 x 10"^ for a random 
fluctuation to be detected as a source ("false-positive"). The first 
stage source lists include 411 and 196 sources for the 2-10 and 
5-10 keV bands, respectively. The source identification numbers 
(ID210 for 2-10 keV and ID510 for 5-10 keV) were assigned at 
this stage. 

For the second stage we used EMLDetect, a standard tool 
from the XMM-Newton S AS softw are. Originally d eveloped for 
ROSAT ([Cmddaceetal]l988i Ha singer et al.|1993| ), the current 



^ The XMM-Newton Users HandBook. |http://xmm.esac.esa.int/| ex- 
ternal/ xmm_user_support/ documentation/uhb/XMM _UHB.html 

^ The correspondence between probability and "number of sigmas" 
only applies for Gaussian probability distributions. Our use should just 
be regarded as a handy mnemonic. 



version of EMLDetect|j includes many improvements in terms 
of PSF models, ability to simultaneously fit data from diff'er- 
ent cameras and pointings, deblend sources, and fit extended 
sourceo 

EMLDetect gives more accurate estimates of net source 
counts, with lowe r systematic and random uncertainties 
( Puccetti et al. '2009). Since EMLDetect was not designed to run 
on a large number of input files, we grouped the observations in 
four sets according to the pointing direction (see Table [T]), in or- 
der to have homogeneous PSFs across a single set. We summed 
together the images from the MOS 1 , M0S2 and PN cameras, ob- 
taining four images which were used as the input to EMLDetect. 
The back ground maps for EMLDet ect were calculated with the 
method of |Cappelluti et al.| ( [2Q07 ). 

We ran EMLDetect using the PWXDetect coordinates as in- 
put and without re-fitting the positions; thus the EMLDetect 
source list is a subset of the PWXDetect one. 

EMLDetect uses a maximum-likelihood estimator to assign 
a significance value to a detection. The EMLDetect source list 
was cut at a likelihood value of L = 4.6. This likelihood obeys 
the law L = -ln(p), where p is the probability of a false-positive 
detection. The threshold L = 4.6 corresponds to a probability of 
1.01 X 10"^. The number of sources detected at the second stage 
is 337 and 135 for the 2-10 and 5-10 keV bands, respectively. 

After visual inspection, we identified two cases of source 
blending, which correspond to ID210s 228 and 276. These 
blends were detected in the 2-10 keV band at 29.1 cr and 22. 3 cr 
by PWXDetect, and with likelihoods L = 565.5 and 340.1 
by EMLDetect, respectively. In both cases, two sources whose 
cores can be visually separated were identified by PWXDetect 
as a single one. The separation between the two components 
is 11.5'' and 18.1'' for ID210s 228 and 276, respectively (see 
also Sect. [7]). Both blends are detected also in the 5-10 keV 
band, as ID510 1109 and 1134, respectively. In the catalogue, 
we do not list the blends; instead we provide counts, rates and 
fluxes estimated with aperture photometry for the four individ- 
ual sources, which in the 2-10 keV catalogue are numbered 501, 
502 (components of 228) and 503, 504 (components of 276). 
These sources have been de-blended also in the 5-10 keV band, 
and their components are numbered with ID510s 1501-1504. 
The main catalogues are defined as the sources detected at the 
second stage, minus the two blends, plus the four blend com- 
ponents, and thus contain 339 and 137 sources in the 2-10 and 
5-10 keV bands, respectively. 

A comparison of the source counts derived by PWXDetect 
and EMLDetect showed a 1-to-l correlation, as expected, but 
with a few outliers. Three of these sources (ID210: 4, 105; 
ID510: 1001) lie at the south-west rim of the field and are proba- 
bly just fluctuations in an area with a strong exposure gradient. A 
few others (ID210: 154, 184, 280, 290, 328, 406; ID510: 1066, 
1089, 1100, 1124, 1146) lie in a crowded area, or very close to 
brighter sources, and they were also flagged as "extended" by 
EMLDetect. Their EMLDetect-derived counts are about one or- 
der of magnitude larger than the PWXDetect estimate, probably 
because EMLDetect fits them as tails of the PSF of a larger and 
brighter source. Aperture photometry showed agreement with 
the counts derived by PWXDetect. Therefore, in the catalogue 
these 6 sources are listed with the PWXDetect counts; count rates 
and fluxes were calculated accordingly. 



^ The XMM SAS reference manual for EMLDetect; 
http://xmm.esac. esa.int/sas/1 1 .0.0/doc/emldetect/node3.html 

^ While the CDFS field containes some extended sources, none of 
them was present in the PWXDetect output. 
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Fig. 4. Colour image of the XMM-CDFS. Red: 
0.4-1 keV; green: 1-2 keV; blue: 2-8 keV. The 
colour scaling is non-linear, and is done first 
by using a quadratic polynomial to sup press the 
background , and then by following the|Lupton| 
|et al.|p004 procedure. 



The flux to count rate conversion factors have be en calcu- 
lated for a power-law spectrum with slope F = 1.7 ("Mainieri" 
et al |2007| ); their values are 1.86 x 10^^ and 0.91 x 10^' counts 



cm~" erg~" for the 2-10 and 5-10 keV bands, respectively. Had 
we used a difl'erent slope (for example, 1.4 or 2.0), the conver- 
sion factors would difl'er by 8-9%. The Galactic column density 
in the CDFS direction is 8 x 10^^ cm"^, so that the bands 2-10 
and 5-10 keV are not significantly aflTected and no correction for 
Galactic absorption is needed. A histogram of the fluxes of the 
XMM-CDFS 2-10 keV catalogue is shown in Fig. (5] The num- 
ber of detected sources drops at fluxes < 10"^^ erg s"^ cm"^, 
with the faintest source at 6.6 x 10"^^ erg s"^ cm"^ giving the 
survey flux limit. For comparison, we also plot the flux distri- 
bution from XI 1 and |Brunner et al.| ( |2008| ). While the Chandra 
4 Ms survey is deeper as expected, the XMM-CDFS contains 
a considerable larger number of sources at medium- faint fluxes 
(2 X 10"^^-2 X 10"^"^ erg s"^ cm"^), most of which are anyway 
detected in the ECDFS. The Lockman hole survey has a distri- 
bution similar to the XMM-CDFS, though a larger number of 
sources (mainly at faint fluxes) is detected in the latter than in 
the former. 



3.1. Catalogue significance 

It is not easy to compute analytically the joint probability of a 
false-positive detection in the two stages, though some hints can 
be obtained using conditioned probability: 



!P(detection) = r{P) P(E\P) 



(1) 



where P(P) is the probability of a PWXDetect false-positive 
detection (6.3 x 10"^), and P(E\P) is the probability of an 
EMLDetect false-positive after the PWXDetect selection. The 
value of ^(detection) can be bounded considering the two ex- 
treme cases of P(E\P) = 1 (i.e., accepting afl PWXDetect 
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Fig.5. 2-10 keV fluxes of the sources detected in the XMM- 
CDFS, compared with other surveys. Solid black histogram: 
main XMM-CDFS catalogue; dashed red histogram: Cha ndra 
catalogue (XI 1); dotted blue histogram: Lockman hole ( ,Brunner| 
let al.|2008| ). 



sources) and P(E\P) = P(P)P(E) (i.e., considering the two de- 
tection stages as independent experiments; P(E) = 1.01 x 10"^): 



P(P)P(E) < !P(detection) < P(P) 



(2) 
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which, in terms of cr, give a threshold significance for the cata- 
logues in the interval Acr-Scr. 

The number of spurious detections can be used to obtain an 
approximation for P{E\P). This number will be first estimated 
in Sect. [6] using simulations, and then refined in Sect. 83 by 
comparing with the Chandra catalogues and by inspecting the 
signal/noise ratios of the sources. Therefore, considering 12 (4) 
spurious sources in the 2-10 (5-10) keV band: 

^detection; 2 - 10) ~ 6.3 x 10"^ x 12/411 ~ 1.8 x 10"^ (3) 
^detection; 5 - 10) ~ 6.3 x 10"^ x 4/196 ~ 1.3 x 10"^ (4) 

which may be quoted as the approximate significance thresh- 
old for the catalogues presented here, and roughly correspond to 
4.11 cr and 4.84(T for the 2-10 and 5-10 keV bands, respectively. 



3.2. Main catalogues 

The XMM-CDFS 2-10 keV and 5-10 keV catalogues include 
339 and 137 sources, respectively, detected in both stages. The 
catalogues are available in electronic form from the Centre de 
Donnees Astronomiques de Strasbourg (CDS), and from the 
XMM-CDFS websitq^ where we will publish updates to the red- 
shifts whenever they become available. The column description 
is for both as follows: 



following 8. 



1. lAUJDENTIFIER — source identifier 
International Astronomical Union conventions; 

2. ID210 (2-10 keV catalogue) or ID510 (5-10 keV catalogue) 
— PWXDetect source number; 
RA — right ascension (degrees); 
DEC — declination (degrees); 
RADEC_ERR — error on position (arcsec; Icr) 
COUNTS — EMLDetect sum of the net source counts from 
M0Sl,M0S2andPN; 

COUNTS _ERR — error on the COUNTS (Icr); 
BKG — (EMLDetect) background counts/arc sec^; 
EXPOSURE — exposure time, averaged over the three cam- 
eras; 

RATE — (EMLDetect) net count rate, summed on the three 
cameras; 
RATE_ERR — error on the RATE (Icr); 



10. 



11. 
12. 
13. 
14. 



FLUX — (EMLDetect) flux (erg s"^ cm'O; 



FLUX_ERR — error on the FLUX (Icr); 
SIGNIFICANCE — (PWXDetect) detection significance 
(no. of sigma); 

15. DETJML — (EMLDetect) detection likelihood; 

16. DET_SCALE — (PWXDetect) wavelet detection scale; 

17. CID — Chandra source number (see Sect.[8]). If the letter E 
appended to the number, then the ID is from [Lehmer et al. 
(2005 ); otherwise, it is from XI 1; 

18. REDSHIFT; 

19. REDSHIFT _REF — reference for choice of redshift; 

20. ID510 (5-10 keV catalogue) or ID210 (2-10 keV catalogue) 
— source number in the other band; 

21. NOTES — mainly about source blends and counterparts. 



3.3. Supplementary catalogues 

The catalogues of the 74 2-10 keV, and 61 5-10 keV sources 
detected only in the first stage are also available in electronic 
form. These supplementary sources are on average detected at 



low significance; many of them are on the borders of the field 
of view; and a few are in crowded fields where EMLDetect has 
trouble separating the diff'erent PSF components. Nonetheless, 
4 of these sources are bright enough that a spectrum could be 
extracted (Comas tri et al., in prep.). 

The columns whose values in the main catalogue were de- 
rived with EMLDetect have been replaced, where possible, 
with analogous quantities derived with PWXDetect. The main 
diff'erences between the output columns of PWXDetect and 
EMLDetect are i) that PWXDetect rescales the MOS counts 
to the PN response before summing them, while EMLDetect 
does not, and ii) that PWXDetect sums the exposure times while 
EMLDetect averages them. For consistency, we renormalized 
the source counts and exposure times to the same scale of 
EMLDetect. 

The column description is as follows: 



1. 



10. 
11. 
12. 
13. 
14. 
15. 
16. 



17. 
18. 
19. 

20. 



lAUJDENTIFIER — source identifier following 

International Astronomical Union conventions; 

ID210 (2-10 keV catalogue) or ID510 (5-10 keV catalogue) 

— PWXDetect source number; 

RA — right ascension (degrees); 

DEC — declination (degrees); 

RADEC_ERR — error on position (arcsec; Icr); 

COUNTS — net source counts from MOSl, M0S2 and PN; 

COUNTS_ERR — error on the COUNTS (Icr); 

BKG — background counts/arcsec^ ; 

SUMMED .EXPOSURE — exposure time, summed on the 

three cameras; 

RATE — net count rate; 

RATE^RR — error on the RATE (Icr); 

FLUX — flux (erg s"^ cm'^); 

FLUX_ERR — error on the FLUX (Icr); 

SIGNIFICANCE — detection significance (no. of sigma); 

DET_SCALE — wavelet detection scale; 

CID — Chandra source number (see Sect. [8]). If the letter E 



appended to the number, then the ID is from Lehmer et al. 

(2005); otherwise, it is from XI 1; 

REDSHIFT; 

REDSHIFT_REF — reference for choice of redshift; 

ID510 (5-10 keV catalogue) or ID210 (2-10 keV catalogue) 

— source number in the other band; 

NOTES — mainly about source blends and counterparts. 



http ://w w w.bo . astro . it/xmmcdf s/deeprime/ 



4. Coverage and number counts 

The survey coverage has been estimated with extensive simu- 
lations, which are described in Sect. [5] In Fig. |6] we show the 
coverage in deg^, defined as the average fraction of sources per 
simulation with a successful match between the output and in- 
put catalogues, times the surveyed area, as function of the flux. 
The confidence intervals for the coverage, estimated according 
to the [Fieller ( |1954 ) formula for the ratio of two average quan- 
tities, are of the order of the line width and thus they are not 
plotted. The coverage is normalized to 0.2466 deg^, i.e. the area 
with exposure >min(r), where T is the exposure column in the 
2-10 keV catalogue. 

In Fig. [6] we also show, for comparison, the coverage for the 
Chandra 4 Ms survey in the 2-10 keV band (XI 1) and that for 
the XMM-Newton Lockman Hole 1 Ms survey for the 5-10 keV 
band ( |Brunner et al.|2008| ). 

The number counts can be derived from the 2-10 and 5- 
10 keV catalogues and the coverage curves derive above. The 
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Fig. 6. Coverage of the XMM-CDFS survey (solid black lines) in the 2-10 keV (left panel) and 5-10 keV band (right panel), defined 
as the fraction of simulated sources which were successfully detected as a function of their output flux, and normalized to the area 
probed by the catalogues. The dotted red line s show, for comparison, the same quantity for Chandra (left panel: 4 Ms survey, XI 1; 
right panel: 1 Ms survey, [Rosati et al.||2002| ). The dashed blue line in the right panel shows the coverage for the XMM-Newton 
Lockman Hole ( Brunner et al.|2008 ). We have not shown in the right panel the CDFS 4 Ms ( Lehmer et al . |20 1 2 ) since formally it is 
in a softer band (4-8 keV). 



cumulative number counts are the sum of the inverse areas in 
which the sources could be detected: 



A^(> 



^)=Ei 



Si>S 



(5) 



With 



where S t and At are the flux and area for a given source - 
this definition, cumulative counts are not binned. 

For the diff'erential counts, we bin the sources according to 
their fluxes, and define 



n{S) 



AS ^ Ai 



(6) 



where the sum is performed on all sources with flux St ^ {S - 
AS 12, S + AS I2\. The Poissonian error on the counts in any bin 
can be computed with the |GehreTs| \ 1 986| ) approximation. 

The number counts are shown in Fig. [7] (upper panels: cumu- 
lative number counts; lower panels: difl'erential). The cumulative 
and diff'erential number counts observed by C handra ([Lehmer 
|etal.|2012| ) and the cumulative model for AGN ('Gini etal.|20"07l ) 
are also shown. The Chandra number counts have been con- 
verted from the 2-8 and 4-8 keV bands to the 2-10 and 5-10 
keV bands, respectively, assuming the same spectrum used for 
the XMM-CDFS catalogue. 

In both bands there is, as expected, a near-perfect agreement 
between the Chandra and XMM-Newton estimates. Some low- 
significance deviations at fluxes around 6 x 10"^^-10"^^ (2-10 
keV) and 2-6 x 10"^^ erg s"^ cm'^ (5-10 keV) might be due 
to the larger field of view of the XMM-CDFS with respect to 
Chandra. 

A cross-correlation analysis between the Chandra ACIS-S3 
detector chi p and the three XMM -Newton cameras in the 2-8 
keV band ( [Tsujimoto et al. 2011| ) showed that Chandra finds 



fluxes which are larger than the XMM-Newton ones by 10-20%. 
The 4 Ms survey however was per formed with the ACIS-I de- 
tector, which was not analyzed by Tsujimoto et al. ( 2011 ). If 
the ACIS-S3 results also hold for the ACIS-I detector, than the 
larger fluxes found by Chandra would be consistent with the mi- 
nor deviations found in the Log N-\.og 5". As to the 5-10 keV 
band, Fig.|7](upper right panel) would suggest that any deviation 
would go in the opposite direction {Chandra fluxes being lower 
than the XMM-A^^w^6>^ ones); however the very hard band was 
not considered by Tsujimoto et al. ( 2011) . 

In both bands the number counts are also in agreement with 
the model. 

The Lockman Hole was also observed with XMM-Newton\ 
Brunner et al.| p008| ) presented a 5-10 keV Log A/^-Log S 
which is in agreement, within poissonian errors, with the XMM- 
CDFS. At brighter fluxes, we compare with the counts from 
the Hellas2XMM (Baldi et al. 2002) and 2XMM (Mateos et al.l 



2008) surveys; all of them are in agreement, within errors, with 
the XMM-CDFS. Poissonian errors are plotted only for the 
CDFS and Hellas2XMM, to avoid cluttering the figure. 

One interesting thing to note is that in the interval 7 x 10~^^- 
1.2 X \{)-^\ the CDFS and the XMM-CDFS 5-10 keV num- 
ber counts stay below all other determinations. This may be a 
low- significance feature of the CDFS field (though formally all 
number counts are consistent with each other within Poissonian 
errors). 

Finally, we show the diff'erential number counts derived 
by [Georgakakis et alT| ( |2008| ) from a combination of deep 
(CDFS and Chandra Deep Field North) and shallow surveys. 
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Georgakakis et al. ( 2Q08| ) report as their best-fit model a broken 
power-law of the form: 



dA^ 
dS 



S^\ S <S 
S^\ S >S 



break 
break 



(7) 



with slopes Ti = -1.56 ± 0.04 and Y2 = -2.52+^-^^, and with 

2-10 



a break flux Log(5'/erg s cm 



-»-2 



keV band (-1.70^^.08^ 



_9 C7+O.O7 
^•^ -0.09' 



) = -13.91^^;^^ for the 



-14.09^QQ5 respectively for 



the 5-10 keV band). The XMM-CDFS data points are con- 
sistent within 1 or 2cr with the broken power-law model. A 
binned x^ fit yields for the XMM-CDFS (quoted errors are 
I0-): Fi = -0.85 ± 1.2, F2 = -2.29 ± 0.36 and break flux 
Log(5/erg s-i cm'^) = -14.46 ± 0.30 for the 2-10 keV band, 
and -0.77 ± 2.6, -2.49 ± 0.48 and -14.47 ± 0.41 for the 5-10 
keV band. T he best-fit paramete r are in agreement, within er- 
rors, with the Georgakakis et al.| p008) ones, with the possible 
exception for the 2-10 keV break flux which difl'ers by 2cr. This 
may be the result of the smaller number of data p oints at fluxes 
fainter than the break in the XMM-CDFS than in Georgakakis 
et al. ( 2008| ), and of the larger error on the faintest data point. 



5. Simulations 

We have performed extensive simulations to address the XMM- 
CDFS properties in terms of survey coverage and source confu- 
sion. For each simulation we built a mock survey, and ran the 
two-stage detection process. 

We have used a simulator developed with the aim to repro- 
duce the background and the PSFs of any given observation as 
closely as possible. Though currently tuned for XMM-Newton, 
the simulator can be adapted to other missions. While more de- 
tails are given in the Appendix, the main features are: 

- the simulator attempts to reproduce the source and back- 
ground distribution for each of the 99 exposures (33 obsids 
times 3 cameras) of the XMM-CDFS independently; 

- the background level is decomposed into its constituents 
(cosmic, particle and residual soft protons), whose spatial 
distributions are reproduced separately; 



- the co smic X-ray background is modeled after [Gilli et al. 
( |2007| ); the particle and proton backgrounds are resampled 
from maps distributed as part of the XMM-Newton calibra- 
tion; 

- the distribution of the XMM-CDFS background counts is re- 
produced by the simulator with an error of ~ 10%; 

- the two-stage source detection is run on each simulation; 

- a sizable number of simulations is run so that average quan- 
tities (number of sources, coverage, number of confused 
sources) can be derived. 

Given the long exposure time (~ 3 Ms) of the XMM-CDFS 
survey, the average number of counts per pixel is much larger 
than that for other XMM-Newton survey projects. Because of 
the difl'erence in the pointing positions among the 33 obsids, the 
background counts have a complex spatial pattern, mainly given 
by the chip gaps and borders occurring in slightly shifted po- 
sitions. In any XMM-A^^w^6>^ observation, the following back- 
ground components are present in addition to the detectable 
sources: 

- cosmic X-ray background (CXB), produced by the many 
sources which are below the survey flux limit; 



- particle background (PB; sometimes also called "non X-ray 
background"), which basically includes the electronic noise 
and a few fluorescent lines produced by the telescope assem- 
bly; 

- residual soft protons (SP): a quiescent background due to soft 
proton clouds; 

- solar wind charge -exchange, negligible at energies > 2 keV 
and therefore not considered here. 

The modelling of each simulated component will be discussed 
in detail in the next Sections. Our treatment of the difl'erent 
background components builds directly on the workby.DeLuca 
|& Molendil ( |2004|), [Snowden et al I (|2004| ) and [Snowden et~al 
( |2008| ); see also ISnowden & KuntzH2012| ). 



5.1. X-ray sources and cosmic background 

For each simulation, an input catalogue of cosmic sources is 
generated from a Log N-Log S obtained as the sum of: i) the 
2-10 k eV Log A^-Log S of AGN, calculated with the [Gilli et aL 
(2007 ) model over its full range of parameters, a nd ii) the 2-10 
keV galaxies number counts computed from the IRanalli et al. 



(2005 ) model. This Log A^-Log S extends from 2.8 x 10"^^ to 
1.0 X 10"^^ erg s"^ cm"^, thus allowing simulation of both the 
detectable sources, and the cosmic X-ray background. It is sam- 
pled on a square area slightly larger than the CDFS, assuming 
uniform distributions in RA and DEC, thus giving a number of 
~ 50000 sources per simulation. 

Most of these sources are too faint to contribute, on average, 
even a single photon to an observation; however when they are 
taken together, the sources below the survey flux contribute to 
the total (cosmic -hparticle-hproton) observed background. 

The flux to count rate conversion factors have been calcu- 
lated for a power-law spectrum with slope F = 1.7, as done for 
the real catalogue (considering F = 1.4 or 2.0 would change the 
fluxes by ~ 8%); in doing this, we have assumed a single aver- 
age spectral slope for all sources, which is acceptable since we 
are interested only in broad-band detections. 



5.2. Particle background 

The PB is mainly due to the instrument electronics and to reflec- 
tion lines from the telescope and its mounting. The spatial dis- 
tribution of the PB has been taken from the Filter Wheel Closed 
(FWC) observations provided by the ESAS CALDB|5 The FWC 
were made with the filter wheel in the closed position, and with- 
out the calibration source 010 

The actual PB level in any observation may be determined by 
considering the corner areas, i.e., the parts of the MOS and PN 
detectors which lie outside of the FOV and are neither exposed 
to the CXB nor to the protons. The data for the corner areas are 
normally filtered out by the XMM-Newton pipeline, but can be 
obtained by reprocessing the "raw data" (i.e., the Observation 
Data Files, ODE). The surface brightness of the corner areas 
needs to be appropriately rescaled to obtain the PB brightness 
inside of the FOV. 

Let us define the total average surface brightnesses which 
are observed inside and outside of the FOV as Ein and Sout, re- 
spectively. The value of both Ss is diff'erent for any obsid and 



^ The Extended Source Analysis Software (ESAS), developed by S. 
Snowden and K. Kuntz, is a part of the standard XMM-Newton SAS 
software and includes several calibration files. 

^ Conversely, the more common CALCLOSED observations do have 
the calibration source visible, in addition to the PB. 
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Fig. 7. Upper panels: cumulative number counts in the 2-10 keV (left) and 5-10 keV (right) bands. Lower panels: differential 
number counts (left: 2-10 ; right: 5-10 keV band) in Euclidean scale (i.e., S^-^dN/dS). For comparison, in all panels we plot the 
Chandra number counts ( jLehmer et al.^2012^ re d circles). In the upp er panels, we also show the cumulative number counts from 
the XMM-A^^w^6>^ survey in the Lockman Hole ( Brunner et al. 12008) blue stars), and the model cumulative cou nts for AGN (|Gilli 
et al.|2007 black dotted line). In the upper right panel, we also plot the number counts from the Hellas2XMM ( Baldi et al.|2002 l 
and 2XMM ( [Mateos et al.|2008| . In the lower panels, we also show the differential number counts from jGeorgakakis et al.| ( |2008] )7 



On the one hand, the FWC only contain PB, therefore: 



^in,FWC - ^PB,FWC 

^out,FWC = k I^in,FWC = k I^PB,FWC 



(8) 



camera, and can also be measured for the FWC. To minimize 
uncertainties due to small number statistics for bright sources, 
the positions corresponding to the 30 brightest sources have been 
masked out during the determination of the above surface bright- 
nesses. This corresponds to removing sources brighter than the 
2 X 10"^"^ erg s"^ cm"^ threshold; at this flux leve l, 43% of the 
background is resolved (using the HEAO-1 value, [Gruber et al.| where k accounts for spatial variations of the PB surface bright- 
[T9991 ). ness. 

On the other hand, for a generic observation, the FOV will 
contain PB, proton and cosmic background components, while 



10 
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Fig. 8. Comparison of the real XMM-CDFS (A) with one simulation (B). The colour scale is the same for both panels and maps the 
2-10 keV intensity. 



the corners will only have PB: 

Xin = SpB + ^SP + ^CXB 
PB 



^out — kHp 



which is also shown in Fig. A.l and Table [21 

We assume the value Scxb = 7.12x 10"^ erg s"^ cm"^ deg"^ 
(9) for the CXB brightness, which results from integrating the AGN 



where XpB, ^sp. and Scxb are surface brightnesses. 

The important assumption here is that the k is the same for 
the FWC and for all the observations, which is equivalent to as- 
suming that the spatial distribution of the PB is the same in all 
observations. Taking k from Eq. ([8]), from Eq. ^ we get 

_ ^in,FWC .^^. 

^out,FWC 

The average fractions of PB with respect to the total of 
CXB-hPB-hprotons are reported in Table |2] The number of PB 
photons to be simulated is XpB xAxE where A is the FOV area 
and E the exposure time. 

5.3. Soft protons 

Clouds of soft protons are often encountered by XMM-Newton, 
especially when the satellite is close to perigee, and cause strong 
background flares. When such a flare is detected, the time inter- 
vals in which it occurs are usually excluded from further anal- 
ysis. However, there is also a quiescent component of soft pro- 
ton background, with much smaller brightness, that may occur 
for the entire length of an observation, and for which the stan- 
dard recipe of identifying and excluding high-background pe- 
riods is not applicable — or which may remain even after the 
flares are excluded. The quiescent soft protons are subject to the 
telescope vignetting, though to a lesser extent than X-ray pho- 
tons (see Fig. 1 7 in Kuntz & Snowdenj2QQ8j . Both .De Luca &| 
|Molendi| ( [2004| ) and jSnowden et al.| ( |2004| ) have considered this 
problem and we refer to these papers for further details. 

Here we estimate the residual soft proton level from Eq. ([9]): 



and galaxy number c ount s from zero to the 2 x 10 ^^ erg s 



cm threshold (Sect. 5.2). The number counts are normalized 



to the HEAO-1 value for the X-ray background ( [Gruber et al. 
1999); had we normalized to any of the higher values derived 
from XMM-New ton data (e.g., jLumb et aL]|2002| or |De Luca| 
& Molendi|2QQ4 ), this number would have been proportionally 
larger. However, the increase would have been compensated by 
a decrease in the soft proton component, with the net efl'ect of 
a slight change in the amount of background vignetting. This 
efl'ect is probably within the uncertainties with which the back- 
ground brightness is measured. 

We now can determine the number of photons contributed 
by the soft protons as Ssp xAxE, where A and E are again the 
FOV area and the exposure time. 

A simulation produced with this recipe is shown in Fig. [8] 
where an image of the real XMM-CDFS is shown with the same 
colour scale for comparison. 



Table 2. Average percentage contribution of the difl'erent back- 
ground components to the overall level. 



Camera 


Band (keV) 


CXB 


Particle 


Proton 


MOSl 


2-10 


1.1+0.4 


81 + 10 


19 + 10 


M0S2 


2-10 


1.3 + 0.4 


80 + 9 


18 + 9 


PN 


2-10 


2.2 + 0.8 


90 + 6 


7.4 + 6.0 


MOSl 


5-10 


0.5 + 0.2 


88 + 8 


12 + 8 


M0S2 


5-10 


0.6 + 0.2 


87 + 8 


12 + 8 


PN 


5-10 


1.5 + 0.5 


95 + 4 


o Q+4.4 
•^•^-3.9 



Notes. All values are per cent. When calculating the brightness levels, 
the areas around the 30 brightest sources have been excluded. Thus, the 



^SP = ^in r ^CXB 



(11) 



cosmic background refers to sources with fluxes < 2 x 10 



erg s 
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Fig. 9. Distribution of the counts per arcsec^ in the simulations 
(average over 5 simulations), compared with the XMM-CDFS, 
for the 2-10 and 5-10 keV bands. The vertical axis shows the 
area covered by each bin. 

5.4. Simulated fields and catalogues 

A comparison between the distribution of counts reproduced by 
the simulator and that of the XMM-CDFS is shown in Fig. 9j 
where the histograms of the counts per arcsec^ are plotted in the 
2-10 and 5-10 keV bands. The counts include both the back- 
ground (mainly the peak of the distribution) and the sources (the 
right tail). In both bands, the agreement seems good, with a hori- 
zontal shift between simulations and the XMM-CDFS of ~ 10% 
around the peaks. 

We ran about 150 simulations in the 2-10 and 5-10 keV 
bands, respectively, to assess the survey properties. For each 
simulation, the two-stage detection process was performed in 
the same way as for the (real) XMM-CDFS main catalogue. 
The output catalogues were cross-correlated with the (input) 
list of simulated sources. For each detected source, its input 
counterpart is defined as the one which minimizes the score S 
( [Cappelluti et al.|2009| ): 



S = 



ARA 

O-RA 



ADEC 

O'DEC 



ARATE\ 



C'"RATE 






(12) 



where A(RA, DEC, RATE) are the diff'erences between the out- 
put and input right ascension, declination, and count rate, respec- 
tively, and the cr are the corresponding errors. 

We show in Fig. [TO] the histogram of the spatial separation 
between the output and input positions. Most sources show a 
displacement of ~ 1.25'', which was already highlighted by 
jCappelluti erSr|p009| ) and |Pineau et al.|p011| ), and therefore it 
seems to be a characteristic ofXMM-Newton. The displacement 
does not occur in any particular direction and occurs even if the 
source detection is run on just one single camera, or on just one 
obsid. The number of matches drops sharply for displacements 

The input and output fluxes of the simulated sources are 
shown in Fig.[TT} only sources with displacements < 6'' are con- 
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Fig. 10. Spatial separation between the input and output posi- 
tions of simulated sources. 

Table 3. Number of spurious sources in the 2-10 keV band sim- 
ulations. 



Exp. CDFS 


Total sim. 


W. count. 


Spurious 


Fract. 


any 335 
> 300 302 


358 + 16 
306 + 15 


328 + 17 
296 + 15 


30 + 5 
10 + 3 


8.4% 

3.3% 



Notes. The column "CDFS" gives the number of sources found in the 
XMM-CDFS which respect the same selection of the simulations (i.e. 
extent< 0.5 and no de-blending; see text). The columns "Total sim.", 
"W.count.", "Spurious" give the average and standard deviation of the 
number of sources per simulation, the number of sources with a coun- 
terpart in the input catalogues, and the number of spurious, respectively. 
The column "Fract." gives the fraction of spurious sources over the total 
number. Exposure thresholds are expressed in ks. 



sidered here. A few outlier per simulations are present, which is 
compatible with the amount found in the XMM-CDFS (Sect. [3]). 
These outliers are mainly due to EMLDetect identifying some 
sources, especially in crowded areas, as extended sources. In the 
compute of the coverage, we identify the outliers has having an 
EMLDetect flux larger than 3 times the PWXDetect flux; the lat- 
ter is then used for the calculation. At bright fluxes (> 10"^"^ erg 
s"^ cm"^), the output values are systematically larger by about 
20% than the input ones. This uncertainty does not afl'ect the flux 
range where most of the sources (82%) are detected. 



6. Number of spurious sources 

We estimate the number of spurious sources in the 2-10 keV 
band by running the detection process on the simulated observa- 
tions. The total number of sources detected in the simulations is 
reported in Table [3] In the 2-10 keV band, the simulations pro- 
duce catalogues which contain on average 358 ± 16 src sim"\ 
i.e. a larger number of sources than the XMM-CDFS, though 
the difl'erence is less than twice the simulations' standard de- 
viation; in around 10% of the simulations a number of sources 
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Fig. 11. Input vs. output fluxes of the simulated sources, from 
100 2-10 keV simulations. Only sources with a match within 
6'' are considered. Some outliers are present, due to EMLDetect 
uncorrectly associating the source with extended emission es- 
pecially in crowded regions, and for which we rather use the 
PWXDetect fluxes; the number of outliers is a few per simula- 
tion. The red line shows the identity relation. At bright fluxes 

(> ^' 

than the input ones 



10 ^^ erg s ^ cm ^), the output values are about 20% larger 



lower than or equal to the XMM-CDFS is detected. The number 
of detected sources is in general dependent on the background 
level; residual variations between the simulations and XMM- 



CDFS background (see Sect. |5.4[ ) may help to explain the dif- 
ference observed in the 2-10 keV band. 

In the Chandra-COSMOS survey, which used th e same de- 
tection software as this paper, [Puccetti et al.| p009| ) found that 
a large number of spurious sources, only present in simulations, 
were characterized by source sizes smaller than the PSF sizes 
— a clear indication of fluctuations detected as sources — and 
removed them. PWXDetect uses the extent column to express 
the source size in "PSF units" (i.e., adimensional); in particular, 
all sources with extent< 0.5 are likely to be spurious. Thus, in 
all simulations we filtered the catalogue to remove source candi- 
dates satisfying this criterion. 

Spurious sources may arise mainly from background fluc- 
tuations, or from an incorrect reproduction of the spatially- 
dependent background level. For instance, the simulations made 
use of spatial maps of the soft proton background, which were 
derived by time filtering of highly flared observations, and may 
still contain small local enhancements originally due to astro- 
nomical sources. Local enhancements are also present in the 
particle background maps, and further investigation would be 
needed to determine wether these are representative of any given 
observation. In general, any background enhancement in the pro- 
ton and particle maps is likely to appear in many simulations, 
since these maps were used as probability distributions with- 
out smoothing on significantly large scales (e.g., larger than a 
few PSF cores). The number of this kind of spurious sources is 
non-trivial to estimate; for this reason, the numbers of spurious 



sources quoted in Table [3] should be regarded as upper limits if 
appHed to the XMM-CDFS; and in Sect. [9] the comparison of 
the Chandra and XMM-Newton catalogues will be used to fur- 
ther refine the estimate of the spurious fraction. 

A check on the position of the candidate spurious sources 
shows that about half of them (i.e. 20 ± 4 in the 2-10 keV band) 
lie in the outer parts of the field, where the exposure is < 300 ks. 
Under the latter threshold, the number of sources detected in the 
simulations with a counterpart in the input catalogues, and the 
number of sources in the XMM-CDFS, match more closely. 

The main reason for the larger fraction of spurious sources at 
low exposure (< 300 ks) is that fluctuations are more important 
at the borders of the FOV than in the inner area — or, in other 
words, that fluctuations are "averaged out" when the exposure 
is large. The spatial distribution of spurious sources has well- 
defined peaks, falling in the low-exposure areas, which can be 
traced back to localized high background features. 

For the rest of this Sect., we define spurious sources as those 
detected in the simulations, but without any counterpart in the 
input catalogues within 6'' from their position; the numbers are 
reported in Table [3] The threshold of 6'' w as chosen after the 
check for Chandra counterparts (Sect. 8.2), in which a tail of 
reliable associations is present up to 6'' especially for sources at 
large off'-axis angles. For the 2-10 keV, the average number of 
spurious sources is 30+5 src sim"^ The number of XMM-CDFS 
sources to which the comparison can be done is 335 (337 sources 
detected at the second stage, without de-blending sources no. 
501-504, and with 2 sources not considered because they have 
extent< 0.5). 

Confused sources may also give rise to spurious sources, in 
the case that a pair of input sources is detected as a single one, 
and that the detected position is distant more than 6'' from both 
the components of the input pair (see Sect. [Tj "double-detection 
pairs"). The number of this kind of spurious sources is ~ 3 per 
simulation, and does not depend on the position in the field. 

A check on the number of spurious sources, based on the 
probability of association with Chandra counterparts and on the 
signal/noise ratio in the XMM-Newton images, is discussed in 
Sect.lO 

The fraction of spurious sources is higher in the 5-10 keV 
band than in the 2-10 keV, probably because localized back- 
ground features contribute with a larger proportion to the to- 
tal counts. Improving this would require a deeper study of 
the XMM-Newton background and an update of the simulator, 
which we defer to a future paper. 



7. Source confusion 

Because of the size of the XMM-CDFS average angular reso- 
lution (8.5'' FWHM, which is the median across the FOV; see 
Sect. [23]) and the faint flux limit of the XMM-CDFS, it is impor- 
tant to understand how the detection process treats close sources. 
Depending on the distance of its components, a pair of sources 
may be mistakenly detected as a single source ("source confu- 
sion"); the ofl'-axis angle dependence of the PSF and local small- 
scale variations of the background also play a role. This may 
also impact the estimated number of spurious sources, if the 
detected position is distant enough from that of the pair com- 
ponents. Here, we use the simulations to estimate the amount 
of source confusion in the XMM-CDFS survey, and compare at 
end with Chandra data. In this Section, we only consider the 2- 
10 keV band; the 5-10 keV range is expected to be less afl'ected 
by source confusion. 
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Fig. 12. Left panel: Number of pairs of sources, as function of the flux ratio between pairs of sources. Sources with fluxes 10 times 
fainter than the local (i.e. position-dependent) flux limits are included. Solid histogram: pairs of input sources. Dotted histogram: 
pairs detected as a single output source. Dashed histogram: pairs where both components have a detection. 

Right panel: Number of pairs of sources, as function of the spatial separation between the two pair components. Line styles as in the 
left panel. Only pairs with flux ratios lower than 3, and whose components are brighter than 6.6 x 10"^^ erg s"^ cm"^ are considered. 



The larger number of input pairs with respect to output ones in the 3 ''-10'' 
centroid and the position of the detection is larger than 6" ; this happens for - 



range is due to pairs where the distance between the 
3 sources per simulations. 



We start our analysis by considering input sources with 
fluxes greater than 10% that of the faintest source detected in 
the real XMM-CDFS, and scaled for lower sensitivity areas 
through an exposure map (i.e., 0.1 x 6.6 x 10"^^/Exp(RA, DEC) 
erg s"^ cm"^, where Exp(RA, DEC) is normalized so that 
max(Exp(RA, DEC)) = 1). This avoids cluttering the results 
with sources too faint to have a significant impact on the detec- 
tion process. For simplicity, we also do not consider triplets of 
sources, nor groups with a larger number of members. With these 
constraints, we identify an average of 61 (108) pairs per simula- 
tion with separation < 10" (< 15") in the input catalogues. 

The number of pairs per simulation, as function of the flux 
ratio between the two components of the pair, is shown in Fig. [12] 
(left panel) as the solid histogram ("input pairs"). Most pairs 
have a flux ratio < 10 between the two components. Pairs with a 
larger difl'erence between the fluxes of the components are con- 
siderably more rare. 



The right panel of Fig. 12 shows the number of pairs as func- 
tion of the separation between the pair components. To be more 
representative of what is observable in practice in the XMM- 
CDFS, we restrict this panel to pairs with a flux ratio < 3, and 
we also drop any pair with a component fainter than the XMM- 
CDFS flux limit (6.6 x 10"^^ erg s"^ cm'^). In the following 
we will refer to this sub- selection as the "detectable pairs". The 
number of pairs with very small separation (< 2") is quite small 
(-0.4 pairs/sim). It grows rapidly at larger separations, stabi- 
lizing at > 8" with a value of 1.7 pairs/sim/arcsec (i.e., every 
arcsec of separation adds 1.7 pairs). 

We consider the following possibilities for the detection of a 
pair of input sources: 

1. detected as two distinct sources; 

2. detected as one single source. 



7.1. Pairs detected as two distinct sources 

We have searched for source pairs where both components 
are individually detected in the output catalogues of simulated 
sources ("double-detections"). 

The number of double-detection pairs, when seen as a func- 
tion of the pair flux ratio (Fig.[T2J left panel, dashed histogram), 
roughly follows the number of input pairs, though it stays lower 
by ~ 1.5 orders of magnitude and has a sharper drop (a factor 
10 vs. a factor 2 for the double-detections and input pairs, re- 
spectively, at flux ratios < 10). Thus, double detections are more 
probable when the flux ratio is small. 

While the number of detectable input pairs rises with the sep- 
aration from to 8" and remains constant afterwards, the num- 
ber of double-detections is consistent with up to 11" and only 
starts to rise afterwards (Fig. 12 right panel). This means that 
input pairs are not recognised as such by the detection process if 
the separation is < 11". 



7.2. Pairs detected as a single source 

Input pairs detected as a single output source are what is usu- 
ally referred to as confused sources. We have searched for them 
by requiring: i) one output source within 6" from the pair cen- 
troid; ii) no other output sources within 6" from the positions 
of the two input sources. The pair centroid is defined as the 
weighted mean of the pair coordinates, using the count rates as 
the weights. The number of single-detection pairs is also shown 
in Fig. 12 (both panels, dotted histograms). By considering only 
pairs with separation < 15 arcsec, their average number is 13 
pairs/simulation. When comparing with the XMM-CDFS (see 
Sect.|8.2|), this number should be taken as a lower limit, since the 
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simulations do not include clustering; in Sect. |8.2[ the number of 
XMM-Newton sources with two or more Chandra counterparts 
is estimated to be 21. 

Almost all pairs with a large flux ratio (> 10) are detected as 
a single source (Fig. 12 left panel), probably because the fainter 
source cannot be distinguished from a fluctuation in the wings 
of the PSF of the brighter source. Conversely, at small flux ra- 
tios the fraction of single-detections drops, and < 1/4 of the in- 
put pairs with flux ratio < 2.5 are single-detected. The reason is 
probably that it is easier to have a double-detection at small flux 
ratios, where the cores of the two PSFs can be better recognised. 
Also, part of the explanation might involve the inclusion of very 
faint sources in the input list. 

Considering the detectable pairs, the number of single- 
detections as function of the pair separation (Fig.[T2J right panel) 
peaks at ~ 9 arcsec, and it always stays lower than the number of 
input pairs, probably because our selection of "detectable pairs" 
is conservative and allows for pairs where one of the components 
is just below the local flux limit. At separations > 9'' the number 
of single-detections starts to decline, and matches the number of 
double-detections at 14''. The total number of single-detections 
with separation < 20'' is 14, which may taken as an estimate of 
the expected number of confused sources in the XMM-CDFS. 

Summarizing, sources closer than 11" cannot be individually 
detected. Sources separated by larger angles may be individually 
detected, the odds for it increasing with the separation, though 
the PSFs may still be blended. Confusion is no longer an issue 
at angles larger than ~ 16". 

8. Chandra and optical identification and redshifts 

In the catalogue tables we include the most probable association 
between our catalogue and the Chandra 4 Ms and ECDFS ones, 
determined with the procedure explained in Sect. |8.1| (blends and 
other possible associations are reported in the NOTES column). 
For the optical identification and for the redshift determi- 
nation, we rely mainly on the work done by XI 1, and report 
the redshift of the optical counterparts from this catalogue. To 
search for more r edshifts and for updat es, we checked spectro - 
scopic redshif ts in Treister et al. ( 2009b} , [Silverman et al. pOlO ), 
[Cooper et al.|p011|), and|Kurk et al.| ( |2012|); and pho tometric 
redshifts in jCardamone et al.| ( |201()i rand |Raflferty et al.| ( [2Qlli) . 



the positional errors), and A is the angular density of Chandra 
sources with flux higher or equal to that of the Chandra mem- 
ber of the association (using the number counts in Luo et al.| 
200^^ multiplied by the combined square positional error (i.e., 
A = ((t| -h (T^) X N(> S)). Seco nd, to obtain the pro bability of 
real association we use Eq. 11 in |Pineau et al.| ( [2QTT] ). 



P(H,. I r) 



1-h 



P(Hc.) 



l-P(Hcp) 



LR(r) 



-i\" 



(14) 



where P(Hcp) is the (unknown) prior probability of a Chandra 
source to be the counterpart of an XMM-Newton source. Given 
the difl'erent characteristics and depths of the XI 1 and ECDFS 
catalogues, we assume P(Hcp) with two constant values: one 
if the source comes from XI 1, another if it comes from the 
ECDFS. 

The (unknown) number of real Chandra counterparts is the 
sum of the probabilities of association, over all associations: 



Aureal =y/^(i^cp I r) 



(15) 



which can be approximated by using the total number of candi- 
date Chandra counterparts N: 



We finally arrive at 



(16) 



(17) 



combining Eqs. 14 



8.1. Theory 

We have cross-correlated the XMM-CDFS 2-10 keV catalogue 
with the Chandra 4 Ms (XI 1) and ECDFS ( .Lehmer et al.|2005| ) 
catalogues. We have used the likelihood-rat i o estim ator with 



a simple prior as described in Pineau et al. ( 201l|p Briefly, 



this method works in two steps. First, for each XMM- Chandra 
source association we define the likelihood-ratio between the 
probability of the association being real vs. it happening by 
chance: 



-r2/2 



LR(r) = 



lA 



(13) 



(Eq. 9 in Pineau et al. [2011") where both r and A are adi- 
mensional. r is the angular distance between both sources di- 
vided by the combined (in quadrature) positional error (i.e. 

r - dl Jcrj^ + cr^ where d is the angular distance and the cr are 



16] and 15 above. This equation can be re- 
solved iteratively for A^reai (starting with A^j-eai = N/2 for in- 
stance) to get P(Hcp) and hence P(Hcp \ r) for each association. 

8.2. Cross-correlation of XMM-Newton and Chandra sources 

We have used the method outlined above to cross-correlate inde- 
pendently our XMM-Newton 2-10 keV and 5-10 keV catalogue 
with the Chandra 4 Ms (using Tables 3 and 6 in X I 1) and the 
ECDFS (using Tables 2 and 6 in Lehmer et al.|2005| ) catalogues. 

Then the cross-correlations were merged; when a source was 
present in both XI 1 and the ECDFS, we only considered the 
XI 1 position. The cross-correlation was performed considering 
all XMM-Chandra associations with r < 5. For each XMM- 
Newton source with more than one possible Chandra counter- 
part, we dropped the association whose probability was < 3% of 
that of the most probable one. We flagged initially as "Good" the 
associations with probabilities of association P(Hcp | r) > 90 per 
cent. 

The restrictions in the cross-correlation catalogue above 
were chosen to provide a manageable size, while likely keep- 
ing all possible associations of interest. This restricted catalogue 
was extensively inspected by eye (taking also into account the 
XMM-Newton and Chandra fluxes) to get our definitive cross- 
correlation catalogue. The visual inspection was carried out by 
superimposing contours of signal/noise ratio (SNR) on top of the 
XMM-Newton and Chandra images, to help locate the peaks of 



^^ This method has been implemented in a plugin to Aladin, and is 
available fromhttp://saada.u-strasbg.fr/docs/fxp/plugin/|. 



^ ^ The 2-8 keV number counts were used if the Chandra source was 
detected in that band, otherwise the 0.5-2 keV number counts were 
used. If the source was only an upper limit in both bands, the upper 
limit in the 2-8 keV band was assumed to be the flux of the source and 
the number counts in that band were used. 
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emission. The SNR contours were also the basis for the anno- 
tations on the source blending present in the NOTES column of 
the catalogue. 

As a result of this inspection, 10 Chandra counterparts to 
2-10 keV XMM-Newton sources had their "Good" flag revoked 
(2 counterparts to 5-10 keV sources), because either the cor- 
responding XMNl-Newton source was spurious, the Chandra 
member of the association had only an upper limit to its 2- 
8 keV flux or the counterpart was from the ECDFS and there 
was a better XI 1 counterpart. Conversely, 13 (5) counterparts 
were considered "Good" despite having probabilities of asso- 
ciation below 90 per cent: this was changed because the po- 
sition of the XMM-Newton source from our source detection 
did not coincide with the peak emission on the corresponding 
XMM-Newton image, or blending with nearby XMM-Newton 
or Chandra sources. In all cases we have kept the positional in- 
formation from PWXDdetect. 

Finally, we obtained 336 "Good" XMM-Chandra associa- 
tions for XMM-Newton sources in the main 2-10 keV cata- 
logue (137 in the main 5-10 keV catalogue), corresponding to 
339 (137) unique XMM-Newton sources. Out of the 315 (130) 
XMM-Newton sources with one or more Chandra counterparts, 
295 (124) had a single counterpart, 19 (5) had two and 1 (1) had 
three. 24 (7) XMM-Newton sources did not have any Chandra 
counterpart. 

8.3. New sources 

To identify candidate real XMM -CDFS sources not previously 
detected by XI 1 or Lehmer et al. (2005) with Chandra (as 
opposed to candidate spurious XMM-Newton detections) we 
adopt the following baseline criterion: a source is considered 
real if there is at least a signal/noise ratio (SNR) contour where 
SNR> 3 in the band where it was detected. 

The SNR was computed, in the Bayesian framework, by 
treating a source's count rate c as a random variable constrained 
to be c > 0. Since the background images have more than 10 
counts over the area where the exposure map % > 10^ s, we use 
a Gaussian approximation to the Poisson statistics, so that for a 
pixel (/, j) the probability that a source is in that pixel with count 
rate c is 



Pijic) ^ n 






(18) 



where the (k, I) indices run over a circular aperture of radius W 
around pixel (/, j), Wki is the value of a PSF centered in (/, j) at 
pixel (k, I) (approximated here by a gaussian of FWHM= 8.5''), 
and Tki and Bj^i are the total and background image values. This 
can be rewritten as 



^ijic) ^ 



exp 



-(c-CQ,ijf 

2crif 



-'vVf-f^(;;^) 



c>0 



-t'r 



(19) 



where erfc(x) = 27i~^ J '^ e'^^dt is the complementary error 
function. Pij(c) is normalized taking into account that c > 
and cojj and crtj are functions of tki, ^kh Tki and Bki around pixel 
(/, 7). Since Pij(c) is symmetric around cqjj for c > 0, by analogy 
with the "pure" Gaussian case, the probability of having detected 
a source at pixel (/, j) is then 



SNR, 



2cojj 

-f' 



Pii(c)dc = 



2-2crfc{co,iJ/criJ^f{2)) 
2 - erfc(co,i,7o-,w ^/(2)) 



(20) 



We have drawn the contours in terms of "Gaussian sigmas" 
taking into account that, for example, SNR=3 corresponds to 
SNR,^ = 0.9973. 

Of the 24 (7) sources in the 2-10 (5-10) keV main catalogue 
without a Chandra counterpart, a SNR > 3 contour could be 
drawn on the 2-10 (5-10) keV images for 1 1 (2) of them, namely 
ID210 5, 85, 176, 189, 207, 224, 280, 348, 381, 402, 407, and 
ID510 1098 and 1150. Of this list, ID210 207 and ID510 1098 
correspond to the same source. Two sources from the supple- 
mentary 2-10 keV catalogue (ID210 186, 410) also satisfy this 
criterion and are thus considered real. The SNR> 3 contour of 
source ID5 10 1 150 includes a blend of ID210 304 (from the sup- 
plementary sample) and 306 (main), but it peaks around ID210 
304; ID210 304 has no Chandra counterpart, while 306 does 
(source 443 in XI 1, whose position coincides with a peak of 
0.5-2 keV emission in XMM-A^^w^6>^). Thus we consider ID210 
304 as the counterpart to ID510 1150, and include 304 among 
the candidate real sources. 

We also considered a source as real, if it has at least a SNR > 
2 contour in the same band of its detection, plus a SNR > 3 
contour in the other band. This happens for one source: ID510 
1 149 (not detected in the 2-10 keV band). 

Thus the total number of candidate "new", XMM-Newton- 
only sources, is 15, divided as follows: 11 main 2-10 keV, 3 
supplementary 2-10 keV, 3 main 5-10 keV (two sources were 
detected in both bands). Conversely, the number of candidate 
spurious sources is 13 in the main 2-10 keV sample and 4 in 
the 5-10 keV. Thus, the spurious fraction in the main 2-10 keV 
catalogue is (24 - ll)/339 ~ 3.8%. 

8.4. Chandra sensitivity iimits and optical/IR counterparts for 
the new sources 

Among the 15 candidate new sources, 2 are outside the 4 Ms 
and ECDFS areas (ID210: 407, 410) and were not observed by 
Chandra. Of the remaining 13, 8 are only covered by the ECDFS 
(ID210/ID510: 5, 186, 189, 207/1098, 224, 348, 381, 402), and 
5 have also data from the 4 Ms survey (ID210/ID510: 85, 176, 
280, 304/1150, 1149; though source 280 fafls in an area with a 
strong exposure gradient). The 4 Ms flux limits are fainter than 
the ECDFS ones by about two order of magnitudes, therefore we 
will focus on the 5 sources falling in the 4 Ms area. The source 
positions are shown in Fig. 14 on an XMM-A^^w^6>^ exposure 
map with the Chandra survey boundaries superimposed. 

Visual inspection of Chandra 4 Ms data showed low- 
significance local enhancements in the positions of the new 
sources. The XMM-Newton fluxes are compared to the Chandra 
sensitivity in the soft and hard bands in Table |4] The Chandra 
limits in 2-8 keV band for sources 85, 280 and 119 inside the 
4 Ms area are a factor 2-3 lower than the flux observed by 
XMM-Newton in the 2-10 keV band, while for source 176 the 
Chandra limit is 15 times lower than the XMM-A^^w^6>^ ob- 
served flux. Only for source 304/1 150 the Chandra limit is above 
the XMM-Newton flux. Therefore, sources 85, 176, 280 and 119 
should, in principle, have been detected by Chandra. 

All the sources inside the Chandra-covered areas are hard 
sources, because no signal is present in XMM-Newton data in 
the 0.5-2 keV band. Also, the hardness ratio^for the 5 sources 
inside the 4 Ms area are > 0.3 (Table HI and consistent with an 






obscured spectrum with Nh > 10 cm ^. The Chandra flux lim 



Calculated following 



Park et al. 



([2006|) as the median of the poste- 
rior porability distribution of the ratio (H - S)/(H + S) where H and S 
are the counts in the 2-10 and 0.5-2 keV band, respectively. 
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Fig. 13. K-band and 3.6yum cutouts of the 15 XMM-CDFS sources not previously detected by Chandra (images from the MUSYC 
and SIMPLE surveys). The XMM-Newton contours from the 2-10 keV band (5-10 keV for sources 1098 and 1149) are super- 
imposed in blue. A green circle with radius 10'' is drawn around the XMM-Newton position. The first, third and fifth rows show 
K-band images; the second, fourth and sixth rows show 3.6yum images. From left to right, and from top to bottom: ID210 5, 85, 176, 
186, 189, 207, 224, 280, 348, 381, 402, 407, 410, and ID510 1149, 1150. K-band data are not available for ID210 407 and 410. 
Each cutout covers an area of 40'' x 42". The K-band images have been smoothed with a Gaussian kernel of size 3 pixel (0.75"). 
The images are shown in histogram-equalized scaling, ranging from light (fainter fluxes) to dark brown (brighter fluxes). 
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410 
ECDFS 407 • 
• dD2 r 



Table 4. XMM-Newton fluxes and Chandra sensitivities for the 
5 candidate new sources falling in the 4 Ms area. 



CDFS 4Ms, 




Fig. 14. Positions and ID210/ID510 numbers of the 15 candidate 
real sources not previously detected by Chandra, superimposed 
on the XMM-Newton 2-10 keV exposure map. The exposure 
map is drawn in logarithmic scale to show also the field border 
where the exposure is very low. The areas previously surveyed 
by Chandra are also shown. 



its were calculated for an unabsorbed spectrum with T = 1.4; for 
obscured sources the limits would be higher depending on the 
column density. Thus it is possible that the diff'erences between 
the Chandra and XMM-Newton eff'ective area curves (Chandra 
having a higher drop at high energies than XMM-Newton) may 
favour the detection by XMM-Newton of very hard sources close 
to the nominal Chandra flux limit in the hard band. Consistently 
source 176, which has the largest ratio between observed XMM- 
Newton flux and Chandra limit, is also the hardest source with 
a HR -0.88 (this HR may be produced by a power-law with 
r = 1.7 and Nh ~ 5 x 10^^ cm"^; source 176 is not detected 
in the 5-10 keV band probably because of its faintness, though 
an excess is still visible in the 5-10 keV image). If the hardness 
of these sources is due to obscuration, then the non-detection by 
Chandra in the 0.5-2 keV band w ould also be justified. N one of 
these sources was included in the | Iwasawa et al.| ( |2012| ) paper 
because their PWXDetect significance is below the threshold set 
in that paper (except for source 407 which has however a low 
number of counts). Further study of these sources may help shed 
light on the obscured AGN population responsible for the miss- 
ing fraction of the cosmic X-ray background, and will be pre- 
sented in a further paper of this series (Ranalli et al., in prep.). 

For the 15 new sources, we have searched for identifications 
in the K and 3 .6//m bands from th e MUSYC ( ,Taylor et al.,2009, ) 
and SIMPLE ( [Damen et al.|201l| surveys. The closest possible 
matches are presented in Table |5} in a few case we list more 
than one possibility due to the crowdness of the field. Finally, 
in Fig. [13] we present K and 3.6yum image cutouts in the with 
XMM-Newton contours sumperimposed. 



ID2 1 0/ID5 1 XMM XMM 



Chandra limits 



flux 



HR 0.5-2 keV 2-8 keV 



85 

176 

280 

304/1150 

1149 



1.2 
8.3 
3.0 
0.21 
2.4 



0.22 
0.88 
0.31 
0.30 
1 



0.080 
0.089 

0.22 
0.073 

0.21 



0.40 
0.53 

1.1 
0.42 

1.2 



Notes. XMM-Newton fluxes are for the 2-10 keV band (except for 
source ID510 1149, which is in the 5-10 keV band). Chandra sensi- 
tivity limits in the 2-8 and 0.5-2 keV bands from the 4 Ms survey. All 
fluxes in units of 10~^^ erg s~^ cm~^. Hardness ratio are from PN data. 



8.5. Number of chance associations 

We have used the simulations detailed in Sect. \5\ to assess how 
many of these XMM-Chandra associations are due to chance: 
i.e., we check the reliability of the match, not that of the individ- 
ual sources. We took the output XMM source lists from five of 
those simulations (for each band) and cross-correlated them with 
the CDFS 4 Ms and the ECDFS catalogues following the pro- 
cedure outlined above (cross-correlation with CDFS 4 Ms and 
ECDFS independently, removal of the duplicated sources, merg- 
ing, subselection out to r < 5, removal of further counterparts 
with low relative probability and selection of associations with 
high probability). We then merged the five resulting catalogues 
and calculated the fraction of simulated XMM sources that had 
one or more counterparts with probability above 90 per cent: out 
of the 1793 (685) sources in those five simulations, 71 (30) had 
one "Good" Chandra counterpart and 6 (3) had two, there were 
none with three or more. These correspond to a total of chance 
associations with one or more Chandra sources of 4.3 ± 0.5 per 
cenj^(4.8^Qg). Therefore, for the 2-10 keV main catalogue of 
339 sources, we would expect between about 13 and 16 chance 
associations of sources, while for the 5-10 keV main catalogue 
of 136 sources there would be between about 5 and 8 spurious 
pairs. 

9. Conclusions 

We have presented the observations, data reduction, catalogues 
and number counts of the XMM-CDFS survey. Currently the 
deepest XMM-A^^w^6>^ survey, the XMM-CDFS observations 
pose several challenges in their reduction. The large number of 
independent exposures (33 obsids times 3 cameras) centred on 
the same field was unprecedented for XMM-Newton. 

- The very large time span of the observations allowed us to 
find an increase by a factor ~ 2 in the instrumental back- 
ground in the years 2008-2010 with respect to the years 
2001-2002, and estimate the variations of the various back- 
ground components. 

- A careful study of the instrumental background of XMM- 
Newton during the XMM-CDFS observations was done to 
produce simulated observations, whose properties were de- 
signed to be as close as possible to the real XMM-CDFS. By 
using the count ratio between in-FOV and out-of-FOV data. 



^^ We have estimated the uncertainties in those fra ctions using a 
Bayesian approach and the binomial distribution ( ,Wall & Jenkins |2003] l 
with the shortest 68 per cent condence interval (S. Andreon, priv. 
comm.). 
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and the Filter Wheel Closed (FWC) observations, the back- 
ground was decomposed in cosmic, particle, and residual 
soft proton. The simulated observations reproduce the spa- 
tial details of the background (chip gaps, a missing MOSl 
CCD, vignetting), and include detected sources distributed 
according to the |Gilli et aT] p007| ) and [Ranalli et al.|pQ05l ) 
model Log A^-Log S. Mock catalogues were produced by 
running the source detection procedure on the simulated ob- 
servations, and were used to calculate the coverage and the 
number of chance XMM-Newton- Chandra associations. 

- We derived the catalogues in the 2-10 and 5-10 keV bands 
with a two-step procedure. First, the PWXDetect wavelet 
software was used to identify candidate sources with a signif- 
icance equivalent to 4cr, and to find their coordinates. Next, 
we used the SAS EMLDetect tool to further check the sig- 
nificance of the sources, and obtain counts, count rates and 
fluxes. The final catalogues contain 339 and 137 sources in 
the 2-10 and 5-10 keV bands respectively, and have a sig- 
nificance equivalent to 4.8cr. 

The faintest sources have fluxes of 6.6 x 10"^^ and 9.5 x 10"^^ 
erg s"^ cm"^, respectively. The flux limits at 50% of the 
sky coverage are 1.8 x 10"^^ and 4.0 x 10"^^ erg s"^ cm"^, 
respectively. The catalogues were cross-correlated with the 
Chandra 4 Ms (XI 1) and ECDFS ( [Lehmer et al.|2005| ), us- 
ing a likelihood-ratio method. Simulations provided an up- 
per limit to the number of spurious sources, whose number 
is better estimated with comparison with Chandra data and 
inspection of the XMM-Newton signal/noise ratio. The spu- 
rious fraction is thus 13/339 ~ 3.8% in the 2-10 keV band, 
and 4/137 ~ 2.9% in the 5-10 keV. 

- Despite the high background level, XMM-Newton was able 
to detect a few extremely hard sources. Further study of these 
objects will help understanding the most obscured AGN pop- 
ulation responsible for the missing fraction cosmic X-ray 
background (Ranalli et al., in prep.). 

- The number counts where derived in both the 2-10 keV 

and 5-10 keV bands, and extend from 6.6 x 10"^^ and 
9.5x10"^^ to 1.1x10"^^ and 6.7x10-^4 ej^g g-i cjn-2j,^gp^^_ 

lively. They are in agreement with previous derivations with 
the Chandra and XMM-Newton telescopes, and in difl'erent 
fields (the CDFS, the Lockman Hole, and the Hellas2XMM 
and 2XMM surveys). The XMM-CDFS currently reaches 
the faintest flux limit in the 5-10 keV band (9.5 x 10"^^ erg 
s"^ cm"^), ~ 3 times fainter than the Lockman Hole (the 
Chandra 4 Ms survey reach es fainter fluxes, tho ugh the band 
is formally softer, 4-8 keV, [Lehmer et al.|2012| ). 

- Finally, the simulations were used to study the source con- 
fusion in the XMM-CDFS relative to the 2-10 keV band. 
An average of 14 confused sources per simulation (i.e. de- 
tected as a single source, but corresponding to two input 
sources separated by < 15'') was found; when comparing 
to the XMM-CDFS, this number should be considered as a 
lower limit because the simulations do not include source 
clustering. The cross-correlation with Chandra yielded 20 
XMM-Newton sources associated with two or three Chandra 
counterparts. 



Appendix A: The XMM-CDFS simulator 

We have developed a general purpose simulator of X-ray CCD 
observations, tuned for the XMM-A^^w^6>^ EPIC camera, but eas- 
ily extensible to other missions. Although its primary use has 
been the validation of the XMM-CDFS catalogue, the simula- 
tor has in fact been written with future missions in mind (from 



eROSITA to concepts such as the Wide Field X-ray Telescope, 
and the project formerly known as XEUS/IXO/Athena). The 
simulator reproduces both point-like sources and the (cosmic, 
instrumental, and soft proton) background. As to sources and 
cosmic background, the only assumptions are that a library of 
position- and energy-dependent PSFs is available, that the ef- 
fective area of the telescope and detector is known, and that an 
exposure map is provided to account for vignetting and eventual 
chip gaps. As to the instrumental background, the levels for the 
particle and soft proton contributions are needed, which in our 
case have been estimated by analysing the existing XMM-CDFS 
observations. 

The main features are: 

- it produces event files, to be analyzed with common X-ray 
data analysis software; 

- it supports an arbitrary number of sources; 

- arbitrary spectra can be assigned to the sources; 

- it calculates the most appropriate PSF for each source, ac- 
cording to position and energy, by interpolating from the PSF 
library. 

The simulator is written in the Perl langu age, making 
use of the Perl Data Languagq^ libraries (PDL; 'Glazebrookj 
|& Economou||1997 ) for numerical computation and FITS in- 
put/output. The Ftools package is used to process the simulated 
event file headers and make it readable by the SAS; the SAS tool 
merge is used to join the source and background event files. 

The computing time for a single-camera, single-obsid simu- 
lation is around 4-5 min on a 2.80 GHz Xeon processor running 
Linux; the time needed to simulate the whole XMM-CDFS on 
a single CPU is therefore around 7.5 hours. The simulations of 
diff'erent cameras/(9Z7^/(i can however be run in parallel. 

The simulator is released under the terms of the General 
Public License (GPL) as published by the Free Software 
Foundation, and is available on the XMM-CDFS websita^ 



A.1. Point-like sources and cosmic X-ray background 

The goal is to produce an event list with the same properties 
(camera, pointing, exposure) of a real observation. This stage 
of the simulator acts basically as a PSF sampler. Given a list 
of coordinates and count rates, for each source it calculates the 
correct PSF according to the off'-axis angle and the spectrum, 
and samples it for a number of photons depending on the count 
rate and exposure time. The PSF library is provided, for XMM- 
Newton, by three calibration files (CCF), one for each of the 
MOSl, M0S2 and PN, which contain normalized images of the 
PSFs at 11 energies (from 0.1 keV to 15 keV) times 6 off'- axis 
angles (from 0' to 15'). Photon times and energies are also as- 
signed to the events, to obtain an output file which looks like a 
real XMM-Newton event file. 

The cosmic X-ray background is simulated by generating 
many faint sources according to the number counts models of 
AGN ( [Gilli et al.|2007] ) and galaxies ( jRanalfi et al.|2005i ). 

The simulator accepts the following input: a list of RA, DEC, 
and count rates of the sources to simulate; a model spectrum 
formatted as an ASCII table; an event file, which is used to get 
the information about the boresight and camera and (optional) 
an exposure map. 

The output is an event list in FITS format, which is readable 
by both PWXDetect and the SAS. 



^"^ Available at http://pdl.perl.org 

^ ^ http :// w w w. bo .astro . it/xmmcdf s/deeprime/index . html 
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Fig. A.l. Intensity of the particle and proton components of the 
quiescent background in the 2-10 keV band for the 33 ohsid^ of 
the XMM-CDFS. The obsids are numbered from 1 to 33 accord- 
ing to observation date; the vertical dotted line marks the separa- 
tion between the 2001-2002 and 2008-2010 observations. Filled 
triangles: PN particle background; empty triangl es: P N residual 
soft protons (after light curve screening; see Sect. |5.3| ); filled cir- 
cles: MOSl particle; empty circles: MOSl protons. The M0S2 
camera has values very similar to the MOS 1 . 



The algorithm works as follows. For each source, the ex- 
pected counts are calculated, and rounded if they are fractional. 
The rounding needs to conserve the sum of the counts from all 
sources, in order to correctly reproduce the level of the cosmic 
X-ray background. The model spectrum is sampled, obtaining 
the event energies, which are binned according to the same sam- 
pling scheme of the PSF. Finally, the PSF is sampled, according 
to the event energies; if needed, PSFs at diff'erent off'- axis angles 
are interpolated. The positions of the sampled events are placed 
at the RA and DEC coordinates specified for the source, with the 
correct position angle. 



A.2. Particle and soft proton background 

The input for the simulation of the background are: the levels 
of the particl e and soft-proton components (as an example, we 
show in Fig. |A.1| the levels for the XMM-CDFS), which can 
be estimated from the in/out of FOV count ratio (see Sects|5.2| 



and |5.3| ); an energy inteval; and an event file used to derive the 
pointing coordinates, the position angle and the exposure time. 

The FWC observations, which consist of event files, are 
filtered for the energy interval resampled with the bootstrap 
method (Efron & Tibshirani,1986j to obtain the simulated parti- 
cle background event file. 

The soft proton are simulated using a series of images, part 
of the ESAS CALDB files. These were originally obtained from 
severely ffared XMM-Newton observations, by filtering out the 
low-background periods and retaining only the ffares. The ap- 



propriate image is sampled, producing a simulated soft proton 
event file. 

The two event files of the background components are finally 
merged with that for the source and cosmic X-ray background. 

Acknowledgements. We thank an anonymous referee whose comments have 
contributed to improve the presentation of this paper. We thank F. Damiani for 
support and help with the PWXDetect software, and K. Kuntz, S. Molendi, and 
S. Sciortino for useful discussions. PR thanks R. Di Luca for very valuable as- 
sistance in recovering data from a problematic disc; and A. Longinotti and all 
the XMM-Newton helpdesk team for their kind, focused and prompt support of 
the XMM-SAS software. 

We acknowledge financial contribution from the agreement ASI-INAF 
1/009/10/0 and from the PRIN-INAF-2011. PR acknowledges a grant from the 
Greek General Secretariat of Research and Technology in the framework of 
the program Support of Postdoctoral Researchers. The XMM-CDFS simula- 
tor was developed by PR as per contract with the Astronomy Department of 
the University of Bologna. FJC acknowledges support by the Spanish ministry 
of Economy and Competitiveness through the grant AYA2010-21490-C02-01. 
WNB acknowledges the NASA ADP grant NNXl 1 AJ59G. 

This research has made use of the Perl Data Language (PDL) which pro- 
vides a high-level numerical functionality for the Perl programming language 
( [Giazebrook & Economou; 1997} http://pdl.per l.org). 



References 

Antonucci, R. 1993, ARA&A, 31, 473 

Baldi, A., Molendi, S., Comastri, A., et al. 2002, ApJ, 564, 190 

Brunner, H., Cappelluti, N., Hasinger, G., et al. 2008, A&A, 479, 283 

Cappelluti, N., Brusa, M., Hasinger, G., et al. 2009, A&A, 497, 635 

Cappelluti, N., Hasinger, G., Brusa, M., et al. 2007, ApJS, 172, 341 

Cardamone, C. N., van Dokkum, R G., Urry, C. M., et al. 2010, ApJS, 189, 270 

Comastri, A., RanaUi, P., Iwasawa, K., et al. 2011, A&A, 526, L9 

Comastri, A., Setti, G., Zamorani, G., & Hasinger, G. 1995, A&A, 296, 1 

Cooper, M. C, Yan, R., Dickinson, M., et al. 2011, ArXiv e-prints 

Cruddace, R. G., Hasinger, G. R., & Schmitt, J. H. 1988, in European 

Southern Observatory Conference and Workshop Proceedings, Vol. 28, 

European Southern Observatory Conference and Workshop Proceedings, ed. 

F Murtagh, A. Heck, & R Benvenuti, 177-182 
Damen, M., Labbe, I., van Dokkum, P. G., et al. 2011, ApJ, 727, 1 
Damiani, R, Maggio, A., Micela, G., & Sciortino, S. 1997, ApJ, 483, 350 
De Luca, A. & Molendi, S. 2004, A&A, 419, 837 
Efron, B. & Tibshirani, R. 1986, Statistical Science, 1, pp. 54 
Fieller, E. C. 1954, Journal of the Royal Statistical Society. Series B 

(Methodological), 16, pp. 175 
Gehrels, N. 1986, ApJ, 303, 336 
Georgakakis, A., Nandra, K., Laird, E. S., Aird, J., & Trichas, M. 2008, MNRAS, 

388, 1205 
Georgantopoulos, I., Comastri, A., VignaH, C, et al. 2013, acceptedby A&A, 

arXiv: 1303.5556 
Giacconi, R. & Zamorani, G. 1987, ApJ, 313, 20 
Giacconi, R., Zirm, A., Wang, J., et al. 2002, ApJS, 139, 369 
GiUi, R., Comastri, A., & Hasinger, G. 2007, A&A, 463, 79 
Glazebrook, K. & Economou, F. 1997, The Perl Journal, 5, 5 
Gruber, D. E., Matteson, J. L., Peterson, L. E., & Jung, G. V. 1999, ApJ, 520, 

124 
Hasinger, G., Burg, R., Giacconi, R., et al. 1993, A&A, 275, 1 
Hopkins, R F, Hernquist, L., Cox, T. J., et al. 2006, ApJS, 163, 1 
Iwasawa, K., Gilli, R., VignaH, C, et al. 2012, A&A, 546, A84 
Kuntz, K. D. & Snowden, S. L. 2008, A&A, 478, 575 
Kurk, J., Cimatti, A., Daddi, E., et al. 2012, ArXiv e-prints 
Lehmer, B. D., Brandt, W. N., Alexander, D. M., et al. 2005, ApJS, 161, 21 
Lehmer, B. D., Xue, Y. Q., Brandt, W. N., et al. 2012, ApJ, 752, 46 
Lumb, D. H., Warwick, R. S., Page, M., & De Luca, A. 2002, A&A, 389, 93 
Luo, B., Bauer, F E., Brandt, W. N., et al. 2008, ApJS, 179, 19 
Lupton, R., Blanton, M. R., Fekete, G., et al. 2004, PASP, 116, 133 
Maccacaro, T., della Ceca, R., Gioia, I. M., et al. 1991, ApJ, 374, 117 
Mainieri, V, Hasinger, G., Cappelluti, N., et al. 2007, ApJS, 172, 368 
Mateos, S., Warwick, R. S., Carrera, F J., et al. 2008, A&A, 492, 51 
Menci, N., Fiore, F, Puccetti, S., & CavaHere, A. 2008, ApJ, 686, 219 
Murdoch, H. S., Crawford, D. R, & Jauncey, D. L. 1973, ApJ, 183, 1 
Page, M. J., Stevens, J. A., Ivison, R. J., & Carrera, F J. 2004, ApJ, 611, L85 
Park, T., Kashyap, V. L., Siemiginowska, A., et al. 2006, ApJ, 652, 610 
Pineau, F.-X., Motch, C, Carrera, R, et al. 2011, A&A, 527, A126 
Puccetti, S., Fiore, F, D'Elia, V, et al. 2006, A&A, 457, 501 
Puccetti, S., VignaH, C, Cappelluti, N., et al. 2009, ApJS, 185, 586 



20 



p. Ranalli et al.: The XMM-CDFS survey: 2-10 and 5-10 keV catalogues 



Rafferty, D. A., Brandt, W. N., Alexander, D. M., et al. 2011, ApJ, 742, 3 

Ranalli, R, Comastri, A., & Setti, G. 2005, A&A, 440, 23 

Rosati, R, Tozzi, R, Giacconi, R., et al. 2002, ApJ, 566, 667 

Schwartz, D. A., Murray, S. S., & Gursky, H. 1976, ApJ, 204, 315 

Setti, G. & Woltjer, L. 1989, A&A, 224, L21 

Silverman, J. D., Mainieri, V., Salvato, M., et al. 2010, ApJS, 191, 124 

Snowden, S. L., Collier, M. R., & Kuntz, K. D. 2004, ApJ, 610, 1182 

Snowden, S. L. & Kuntz, K. D. 2012, Cookbook for analysis procedures for 

XMM -Newton EPIC MOS observations of extended sources and the diffuse 

background, version 5, Tech. rep. 
Snowden, S. L., Mushotzky, R. R, Kuntz, K. D., & Davis, D. S. 2008, A&A, 

478, 615 
Taylor, E. N., Franx, M., van Dokkum, R G., et al. 2009, ApJS, 183, 295 
Treister, E. & Urry, C. M. 2005, ApJ, 630, 115 
Treister, E., Urry, C. M., & Virani, S. 2009a, ApJ, 696, 110 
Treister, E., Virani, S., Gawiser, E., et al. 2009b, ApJ, 693, 1713 
Tsujimoto, M., Guainazzi, M., Plucinsky, R R, et al. 2011, A&A, 525, A25 
Urry, C. M. & Padovani, R 1995, PASP, 107, 803 
Wall, J. & Jenkins, C. 2003, Practical Statistics for Astronomers, Cambridge 

Observing Handbooks for Research Astronomers (Cambridge University 

Press) 
Xue, Y. Q., Luo, B., Brandt, W. N., et al. 2011, ApJS, 195, 10 



21 



p. Ranalli et al.: The XMM-CDFS survey: 2-10 and 5-10 keV catalogues 
Table 1. XMM-CDFS observation log. 



OBSJD 


DATE 


RA 


DEC 


R\ 


MOSl EXP 


M0S2 EXP 


PNEXP 


BORES.GR. 


0108060401 


2001-07-27 


53.0812 


-27.7929 


58.74 


^24.6 


^24.6 


« 18.8 


A 


0108060501 


2001-07-27 


53.0919 


-27.7981 


58.75 


43.6 


44.6 


35.8 


A 


0108060601 


2002-01-13 


53.1457 


-27.8252 


238.71 


^52.6 


^52.6 


43.6 


B 


0108060701 


2002-01-14 


53.1406 


-27.8218 


238.72 


78.5 


78.7 


69.4 


B 


0108061801 


2002-01-16 


53.1450 


-27.8195 


238.72 


55.1 


54.9 


54.4 


B 


0108061901 


2002-01-17 


53.1456 


-27.8133 


238.71 


42.9 


42.8 


39.1 


B 


0108062101 


2002-01-20 


53.1511 


-27.8168 


238.72 


44.4 


44.3 


42.6 


B 


0108062301 


2002-01-23 


53.1462 


-27.8147 


238.71 


83.7 


82.5 


72.3 


B 


0555780101 


2008-07-05 


53.1476 


-27.7368 


61.73 


106.6 


109.3 


88.8 


C 


0555780201 


2008-07-07 


53.1488 


-27.7457 


61.72 


116.7 


118.7 


^91.5 


C 


0555780301 


2008-07-09 


53.1387 


-27.7448 


61.73 


n09.0 


108.6 


96.7 


C 


0555780401 


2008-07-11 


53.1385 


-27.7371 


61.74 


89.9 


91.5 


77.7 


C 


0555780501 


2009-01-06 


53.1325 


-27.8371 


241.73 


97.9 


99.2 


90.7 


D 


0555780601 


2009-01-10 


53.1343 


-27.8439 


241.71 


104.9 


108.8 


79.3 


D 


0555780701 


2009-01-12 


53.1322 


-27.8545 


241.72 


101.5 


102.5 


95.7 


D 


0555780801 


2009-01-16 


53.1233 


-27.8359 


241.71 


88.4 


89.1 


77.0 


D 


0555780901 


2009-01-18 


53.1230 


-27.8445 


241.74 


96.0 


96.3 


80.6 


D 


0555781001 


2009-01-22 


53.1234 


-27.8528 


241.73 


106.4 


108.2 


101.3 


D 


0555782301 


2009-01-24 


53.1235 


-27.8531 


241.74 


106.5 


107.1 


96.4 


D 


0604960101 


2009-07-27 


53.1486 


-27.7447 


62.73 


101.5 


105.1 


86.8 


C 


0604960201 


2009-07-17 


53.1404 


-27.7448 


62.73 


104.7 


104.4 


^76.9 


C 


0604960301 


2009-07-05 


53.1496 


-27.7550 


62.72 


noi.4 


^ 102.3 


^89.0 


C 


0604960401 


2009-07-29 


53.1393 


-27.7544 


62.73 


119.2 


120.2 


112.5 


C 


0604960501 


2010-01-18 


53.1307 


-27.8384 


245.73 


46.3 


46.5 


45.0 


D 


0604960601 


2010-01-26 


53.1312 


-27.8463 


245.73 


^ 106.9 


^ 106.8 


95.6 


D 


0604960701 


2010-01-12 


53.1224 


-27.8388 


245.73 


^91.7 


^92.6 


^51.1 


D 


0604960801 


2010-02-05 


53.1315 


-27.8569 


245.73 


90.3 


90.3 


82.3 


D 


0604960901 


2010-02-11 


53.1221 


-27.8573 


245.72 


^84.7 


«84.8 


75.4 


D 


0604961001 


2010-02-13 


53.1211 


-27.8476 


245.72 


91.4 


93.7 


86.1 


D 


0604961101 


2010-01-04 


53.1316 


-27.8301 


245.72 


104.2 


105.9 


100.6 


D 


0604961201 


2010-01-08 


53.1197 


-27.8294 


245.71 


114.0 


115.0 


^96.9 


D 


0604961301 


2010-01-19 


53.1307 


-27.8379 


245.71 


M1.2 


M1.9 


^5.8 


D 


0604961801 


2010-02-17 


53.1213 


-27.8474 


245.73 


94.9 


95.5 


88.8 


D 



Notes. The columns show: the obsid number; the date of observation start; the pointing coordinates (RA and DEC in deg and in J2000 coordinates; 
position angle in deg); the exposure times in ks for the MOSl, M0S2 and PN cameras after background filtering; and the boresight group (see 
Fig. [3]). ^: highly flared observations for which the 3cr-clipping failed (see text). 
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Table 5. K-band and 3.6yum counterparts to new XMM-CDFS sources. 



ID210/ 


Location 


RAnir 


DECnir 


dist 


ID 


K 


ID 


3.6yum 


redshift notes 


IDS 10 




J2000 


J2000 


arcsec 


MUSYC 


AB 


SIMPLE 


AB 




407 


out 


53.25516 


-27.54088 


1.46 


- 


- 


54179 


18.58 


1.13 ^'^ 


410 


out 


53.22289 


-27.51963 


0.83 


- 


- 


56022 


21.33 


_ c 


5 


ECDFS 


53.09970 


-28.02297 


1.92 


1613 


20.56 


11010 


20.31 


0.68 


186 


ECDFS 


53.89292 


-27.80909 


1.09 


8374 


19.71 


30699 


19.67 


- 


189 


ECDFS 


52.90288 


-27.80730 


1.60 


8413 


21.35 


30653 


20.02 


1.43 


207/1098 


ECDFS 


52.88460 


-27.78947 


1.81 


8998 


20.94 


32579 


20.91 


0.62 


224 


ECDFS 


53.32407 


-27.77321 


1.86 


9603 


20.51 


34104 


19.80 


1.38 ^ 


348 


ECDFS 


53.36649 


-27.64710 


3.31 


- 


- 


46367 


23.17 


_ e 






53.36634 


-27.64582 


2.71 


- 


- 


46310 


21.71 


_ e 






53.36575 


-27.64786 


6.89 


13769 


21.81 


46252 


22.63 


- 


381 


ECDFS 


53.03626 


-27.59340 


4.2 


- 


- 


50549 


21.56 


_ e 






53.03574 


-27.59456 


2.40 


- 


- 


50687 


23.64 


_ e 






53.03471 


-27.59362 


2.32 


- 


- 


50766 


23.69 


_ e 


402 


ECDFS 


53.18346 


-27.56021 


3.14 


15482 


20.51 


52934 


19.79 


f 






53.18321 


-27.56246 


7.74 


16344 


18.95 


52602 


19.53 


- 


85 


4 Ms 


53.07197 


-27.90447 


1.93 


4907 


21.37 


21475 


20.53 


1.1125 


176 


4 Ms 


53.25321 


-27.81607 


0.12 


8087 


21.49 


30251 


21.17 


1.75 


280 


4 Ms 


53.26146 


-27.72112 


2.64 


10952 


15.63 


37523 


16.42 


_ d 


1149 


4 Ms 


52.99465 


-27.70056 


3.2 


- 


- 


41386 


23.93 


_ e 






52.99755 


-27.70113 


7.4 


12085 


21.23 


41190 


20.53 


- 


304/1150 


4 Ms 


53.13756 


-27.70008 


2.49 


12132 


21.54 


41087 


20.57 


1.2357 



Notes. The columns show: ID210 (ID510 if the number is > 1000); if the source falls inside or outside the 4 Ms and ECDFS surveys; J2000 
positions from the K band, unless in case of a missing K counterpart, in which case we list the IRAC 3.6//m position; distance in arcsec between 
the XMM-Newton and NIR position; ID number and K magnitude from the MUSYC survey; ID number and 3.6yt/m magnitude from the SIMPLE 
survey; redshift; notes. Magnitudes are in the AB scale. 

(«) No K-band imaging. ^^^ Bright XMM-Newton source. ^''^ Possible extended NIR source, only closest XMM-A^^w^6>;i/NIR match Hsted. ^'^^ K-band 
empty field. ^^^ Another source is present at 7.7''. ^^^ Bright galaxy. 
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